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Abstract 

In this letter, we make a prima facie case that there could be distinct advan- 
tages to exploiting a new class of finite flux equilibrium solutions of the Quantum 
Boltzmann equation in semiconductor lasers. 

1 Introduction 

At first sight, it may very well seem that the two subjects linked in the title have 
little in common. What do semiconductor lasers have to do with behavior normally 
associated with fully developed hydrodynamic turbulence? In order to make the 
connection, we begin by reviewing the salient features of semiconductor lasers. In 
many ways, they are like two level lasers in that the coherent light output is asso- 
ciated with the in phase transitions of an electron from a higher to lower energy 
state. In semiconductors, the lower energy state is the valence band from which sea 
electrons are removed leaving behind positively charged holes. The higher energy 
state is the conduction band. The quantum of energy released corresponds to an 
excited electron in the conduction band combining with a hole in the lower band 
below the bandgap. Bandgaps, or forbidden energy zones are features of the energy 
spectrum of an electron in periodic potentials introduced in this case by the periodic 
nature of the semiconductor lattice. 

However, there are two important ways in which the semiconductor laser differs 
from and is more complicated than the traditional two-level laser model. First, there 
is a continuum of bandgaps parameterized by the electron momentum k and the 
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laser output is a weighted sum of contributions from polarizations corresponding 
to electron-hole pairs at each momentum value. In this feature, the semiconduc- 
tor laser resembles an inhomogeneosly broadened two level laser. Second, electrons 
and holes interact with each other via Coulomb forces. Although this interaction is 
screened by the presence of many electrons and holes, it is nonetheless sufficiently 
strong to lead to a nonlinear coupling between electrons and holes at different mo- 
menta. The net effect of these collisions is a redistribution of carriers (the common 
name for both electrons and holes) across the momentum spectrum. In fact it is the 
fastest (ps 100 fs.) process (for electric field pulses of duration greater than picosec- 
onds) and because of this, the gas of carriers essentially relaxes to a distribution 
corresponding to an equilibrium of this collision process. This equilibrium state is 
commonly taken to be that of thermodynamic equilibrium for fermion gases, the 
Fermi-Dirac distribution characterized by two parameters, the chemical potential \i 
and temperature T, slightly modified by the presence of broadband pumping and 
damping. 

But the Fermi-Dirac distribution is not the only equilibrium of the collision 
process. There are other stationary solutions, called finite flux equilibria, for which 
there is a finite and constant flux of carriers and energy across a given spectral 
window. The Fermi-Dirac solution has zero flux of both quantities. It is the aim of 
this letter to suggest that these finite flux equilibria are more relevant to situations 
in which energy and carriers are added in one region of the spectrum, redistributed 
via collision processes to another region where they are absorbed. Moreover, it 
may be advantageous to pump the laser in this way because such a strategy may 
partially overcome the deleterious effects of Pauli blocking. The Pauli exclusion 
principle means that two electrons with the same energy and spin cannot occupy the 
same state at a given momentum. This leads to inefficiency because the pumping 
is effectively multiplied by a factor (1 — n s (k)),s = e, h for electrons and holes 
respectively, denoting the probability of not finding electron (hole) in a certain k 
(used to denote both momentum and spin) state. But, near the momentum value 
corresponding to the lasing frequency u>l, n s (k) is large (n e (k) + n/j(k) must exceed 
unity) and Pauli blocking significant. Therefore, pumping the laser in a window 
about loq > hi such a way that one balances the savings gained by lessening 
the Pauli blocking (because the carriers density n s (k) decreases with k = |k|) with 
the extra input energy required (because k is larger), and then using the finite 
flux solution to transport carriers (and energy) back to lasing frequency, seems 
an option worth considering. The aim of this letter is to demonstrate, using the 
simplest possible model, that this alternative is viable. More detailed results using 
more sophisticated (but far more complicated) models will be given later. 

These finite flux equilibria are the analogies of the Kolmogorov spectra associ- 
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ated with fully developed, high Reynolds number hydrodynamic turbulence and the 
wave turbulence of surface gravity waves on the sea. In the former context, energy 
is essentially added at large scales (by stirring or some instability mechanism), is 
dissipated at small (Kolmogorov and smaller) scales of the order of less than the 
inverse three quarter power of the Reynolds number. It cascades via nonlinear inter- 
actions from the large scales to the small scales through a window of transparency 
(the inertial range in which neither forcing nor damping is important) by the con- 
stant energy flux Kolmogorov solution. Indeed, for hydrodynamic turbulence, the 
analogue to the Fermi-Dirac distribution, the Rayleigh-Jeans spectrum of equipati- 
tions, is irrelevant altogether. The weak turbulence of surface gravity waves is the 
classical analogue of the case of weakly interacting fermions. The mechanism for 
energy and carrier density (particle number) transfer is "energy" and "momentum" 
conserving binary collisions satisfying the "four wave resonance" conditions 

k + ki = k 2 + k 3 , w(k) +w(ki) = w(k 2 ) +w(k 3 ). (1) 

In the semiconductor context, hu>(k) = hui ga _ p + e e (k) + e^(k) (which can be well 
approximated by a+fik 2 ) where haj ga _ p = e gap corresponds to the minimum bandgap 
and e e (k), e/ l (k) are electron and hole energies. In each case, there is also a simple 
relation E(k) = um(k) between the spectral energy density E(k) and carrier (par- 
ticle number) density n(k). As a consequence of conservation of both energy and 
carriers, it can be argued (schematically shown in Figure 1 and described in its cap- 
tion), that the flux energy (and some carriers) from intermediate momentum scales 
(around ko say) at which it is injected, to higher momenta (where it is converted 
into heat) must be accompanied by the flux of carriers and some energy from fco to 
lower momenta at which it will be absorbed by the laser. It is the latter solution 
that we plan to exploit. 

2 Model 

We present the results of a numerical simulation of a greatly simplified model of 
semiconductor lasing in which we use parameter values which are realistic but make 
fairly severe approximations in which we (a) assume that the densities of electrons 
and holes are the same (even though their masses differ considerably) (b) ignore 
carrier recombination losses and (c) model the collision integral by a differential 
approximation Q , , |J in which the principal contributions to wavevector quartets 
satisfying (|T|) are assumed to come from nearby neighbors. Despite the brutality 
of the approximations, the results we obtain are qualitatively similar to what we 
obtain using more sophisticated and complicated descriptions. 
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The semiconductor Maxwell-Bloch equations are (|]-(5|, 
de ■ n f 

di = l 2e~ J Wkd ~ lEe ' ^ i 

^ = (iQ-iuj k -'yp)p k - l -^-(2n k -l)e, (3) 
dn k ( dn k \ i 



dt 



A(l - n k ) - -y k n k + i-^f] - (p, k p k e* - fi k p* k e) . (4) 

V 01 / collision zn 



Here e(t) and p k {t) are the electric field and polarization at momentum k envelopes 
of the carrier wave exp (—i£lt + iKz) where is the cavity frequency (we assume 
single mode operation only) and n(k) is the carrier density for electrons and holes. 
The constants je, 7p model electric field and homogeneous broadening losses, eo is 
dielectric constant, \i k is the weighting accorded to different k momentum (modeled 
by = Atfc=o/(l + e fc/ e gap)), Afc and 7& represent carrier pumping and damping. In 
(||), the collision term is all important and is given by 

d f 

—n k = 4ir \T kklk2ks \ 2 (n k2 n k3 (l - n kl - n k ) + n k n kl (n k2 + n ka - 1)) 

xd(k + ki - k 2 - k 3 )(5(o; k + w kl - w k2 - w k3 )dkidk 2 dk 3 , (5) 

where T kklk2k3 is the coupling coefficient measuring mutual electron and hole in- 
teractions. We make the weak assumption that all fields are isotropic and make a 
convenient transformation from k (= |k|) to u via the dispersion relation uj = to(k) 
defining the carrier density by J N^duj = j n(k)dk or N w = A7rk 2 dk/(du>)n(k). 
Then, in the differential approximation, (|5|) can be written as both, 

dNu 8 2 K 3uN u 8 8K. 

and -^ = -7rl K - UJ 7T7^ ( 6 ) 



dt duj 2 dt du dtu 



with 



K = 



(nt ( n j) + n% Qnn u )"'\ , ()' = n u = n(k(w)) 



where s is the number computed from the dispersion relation, the dependence of 
2fcfcifc 2 fc3 on k and dimensions (s is of the order of 7 for semiconductors.) The 
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conservation forms of the equations for N w and E u = ojN^ allow us identify Q = 



(positive if carriers flow from high to low momenta) and P = K — uJ^jjj (positive 
if energy flows from low to high momenta) as the fluxes of carriers and energy 
respectively. Moreover, the equilibrium solutions are now all transparent. The 
general stationary solution of @ is the integral of K = Qlo + P which contains four 
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parameters, two (chemical potential and temperature) associated with the fact that 
K is a second derivative, and two constant fluxes Q and P of carriers and energy. 
The Fermi-Dirac solution = (exp (Auj + B) + 1) , the solution of K = 0, has 
zero flux. We will now solve (§), (||) and @ after angle averaging (f|) and replacing 
Ank 2 §r- ( tt^ I by The value of the constant / is chosen to ensure that 

solutions of (||) relax in a time of 100 fs. 

3 Results 

We show the results in figures 2, 3, and 4. First, to test accuracy, we show, in Figure 
2, the relaxation of @ to a pure Fermi-Dirac spectrum in the window lol = 1 < u < 
u>o = 2. The boundary conditions correspond to P = Q = at both ends. We then 
modify boundary conditions to read Q = Qq > and P = at both ends. Next, 
in Figure 3 and 4, we show the results of two experiments in which we compare the 
efficiencies of two experiments in which we arrange to (i) pump broadly so that the 
effective carrier distribution equilibrium has zero flux and (ii) pump carriers and 
energy into a narrow band of frequencies about loq and simulate this by specifying 
carrier and energy flux rates Ql and Pl = —ujlQl (Pl chosen so that the energy 
absorbed by the laser is consistent with the number of carriers absorbed there) at 
the boundary u = ujq. uj = ojl is the frequency at which the system lases. 

In both cases, the rate of addition of carriers and energy is (approximately) 
the same. The results support the idea that it is worth exploring the exploitation 
of the finite flux equilibrium. The carrier density of the equilibrium solutions at 
loq is small thus making pumping more efficient there. The output of the laser is 
greater by a factor of 10. While we do not claim that, when all effects are taken 
account of, this advantage will necessary remain, we do suggest that the strategy of 
using finite flux equilibrium solutions of the Quantum Boltzmann equation is worth 
further exploration. 
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Figure Captions. 



Figure 1 

Carriers and energy are added at ojq at rates Qo and loqQq. Energy and some 
carriers are dissipated at ur > ujq (an idealization) and carriers and some 
energy are absorbed by the laser at ul- (The carriers number will build until 
the laser switches on.) A little calculation shows Ql = Qoi^R — ^ J o)/{ UJ R — 
ul), Qr = Qo(u L - loq)/(lo r - u L ), P R = Q uj r (u;o - lo l )/{ur - u L ), P L = 
^L<5o(^o — ur)/{ojr — ul). Finite flux stationary solutions are realized in the 
windows (wl,wo) and (wo,wr) although in practice there will be some losses 
through both these regions. 

Figure 2. 

To test accuracy we take some initial distribution function (thin line) and plot 
its time evolution as described by (JgJ) with boundary conditions P = Q = 
at both ends. The distribution function relaxes to Fermi-Dirac state (thick 
line). Several intermediate states is shown by long-dashed and short-dashed 
lines (Figure 2a). We then modify boundary conditions to P = 0,Q = Qo > 
at both ends. Then initial distribution function (thick line) relaxes to finite-Q- 
equilibria as shown by long-dashed line. We then change boundary conditions 
to P = 0, Q = Qi > Qo at both ends, so that distribution function is shown 
by short-dashed line. Increasing Q at boundaries even further, so that P = 
0; Q = Qi > Qi a t both ends, the distribution function is given by dotted line 
(Figure 2b). 

Figure 3. 

We now solve (§-|]) with the collision term given by (^). We pump broadly, 
so that the effective carrier distribution has zero flux. The initial distribution 
function (thin line) builds up because of a global pumping (dashed lines), until 
the laser switches on. The final (steady) distribution function is shown by thick 
solid line (Figure 3a). The output power (in arbitrary units) as a function of 
time (measured in relaxation times ~ lOOfs) is also shown (Figure 3b). 

Figure 4. 

We pump in the narrow region around u>o ~ 200meV and we model this by 
specifying carrier and energy flux rates Ql and Pl = —ujlQl- The initial 
distribution function (thin line) builds up because of influx of particles and 
energy from right boundary (dashed lines), until the laser switches on. The 
final (steady) distribution function is shown by thick solid line and corresponds 
to a flux of particles and energy from right boundary (where we add particles 
and energy) to the left boundary, where the system lases (Figure 4a). The 
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output power as a function of time is also shown (Figure 4b). 
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Figure 1 
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Figure 2 (a) 

Relaxation to Fermi -Dirac distribution 
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Figure 3 (a) 
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Figure 4 (a) 
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Figure 4 (b) 
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